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The effects of impurity and spin-orbit scattering potentials can strongly 
affect the Zeeman response of a c?-wave superconductor. Here, both the phase 
diagram and the quasiparticle density of states are calculated within the Born 
approximation and it is found that the spin-orbit interaction influences in a 
qualitatively different way the Zeeman response of d-w&ve and s-wave super- 
conductors. 
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I. INTRODUCTION 

The layered structure of cuprates makes these materials good candidates for observing 
Zeeman response to a magnetic field H directed parallel to the Cu-0 planes |],@]. More- 
over, the d x 2_ y 2 symmetry of the order parameter (hereafter <i-wave) leads in principle to 
substantial differences with respect to the Zeeman response of isotropic s-wave supercon- 
ductors. For example, at zero temperature, the tunneling conductance a s (0) of a (i-wave 
superconductor-insulator-metal junction is nonzero for finite voltages V provided H ^ 
Hill, in sharp contrast to ordinary isotropic s-wave junctions for which a s (0) is zero for 
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V < A/e, where A is the energy gap and e is the electron charge ||. On the other hand, 
the phase diagrams of pure s-wave and ci-wave superconductors in the presence of a Zee- 
man magnetic field have similar qualitative behaviours. For example, for both symmetries 
of the order parameter, a first order phase transition to the normal state is found at low 
temperatures and for sufficiently strong magnetic fields However there are quantitative 
differences. For example, at T = 0, the critical field is HbH c / Ao = l/y/2 for s-wave 
and fj,BH c /Ao ~ 0.56 for <i-wave HUH, where Ao is the zero temperature order parameter 
without magnetic field. 

So far, systematic theoretical studies of the Zeeman response of anisotropic supercon- 
ductors have been focused on the clean limit of <i-wave BCS formulation. A more realistic 
situation would require the inclusion of impurity effects, since these are known to have im- 
portant effects on both thermodynamic and spectral quantities ||. Moreover, in addition 
to the disorder potential, the quasiparticles are also spin-orbit coupled to the impurities, 
so that the Zeeman response is affected by spin-mixing processes. An additional source for 
spin-orbit effects could be provided by the electric fields in the vicinity of the conducting 
Cu-0 layers and the charge reservoirs interfaces. 

Although, in the past years, the effect of spin-orbit coupling has been largely studied 
for isotropic s-wave superconductors the corresponding situation for d-wave supercon- 
ductors (or other anisotropic symmetries) is still unknown. However, it is expected that 
the spin-orbit effects on the Zeeman response of d-wave superconductors differ from those 
of s-wave superconductors in a qualitative way. In fact, already at zero magnetic field, the 
spin-orbit scattering is pair-breaking and reduces both the critical temperature T c and the 
order parameter ||. As a consequence, for H ^ 0, the pair-breaking effects of both the 
external magnetic field and the spin-orbit coupling add together. This situation must be 
contrasted with the s-wave case, where the spin-orbit potential is not pair-breaking and 
competes with the Zeeman response reducing the pair-breaking effect of the magnetic field 

§■ 

In this paper, the effects of both impurity and spin-orbit scattering potentials are stud- 
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ied within a self-consistent Born approximation for ci-wave superconductors. Both ther- 
modynamic and spectral properties are investigated and compared with those of s-wave 
superconductors. 

II. BORN APPROXIMATION 

Let us consider a two-dimensional system with electrons (holes) moving in the x-y plane 
under the effect of an external magnetic field H directed along the plane. In this situation, 
the coupling of the orbital motion of the charge carriers to the magnetic field is vanishingly 
small. In the following, no particular pairing mechanism is assumed and the condensate will 
be described within the BCS formalism. In this framework, the hamiltonian is: 

H = ^e(k)c kQ c ka - lJ2 ac ka c ka ~ A ( k )( c kt c -k| + C-k|C kT ), (1) 

k,a ka k 

where / = hbH and hb is the Bohr magneton. For a d x 2_ y 2 symmetry of the gap, A(k) is 
parametrized as follows: 

A(k) = Acos(20), (2) 

where <fi is the polar angle in the k x -k y plane. In Eq.(|l|), j and j refer to the spin direction 
along and opposite to the direction of H and it is assumed that H is directed along the x 
direction, so that H = i/x. 

The interaction hamiltonian describing the coupling to the impurities located randomly 
at Rj is given below: 

H = Vi m p e ^ ^ l c\ La ci L ' a 

kk'i a 

+ iw £ E e- l(k - k,) - R *([k x k'] ■ tr^ce, (3) 

F kk,i a,P 

where v imp and v so refer to the non-magnetic and spin-orbit coupling to the impurities, 
respectively (kp is the Fermi momentum). From the Elliott relation 0, the impurity and 
spin-orbit potentials are roughly given by v so ~ Agvi mp , where Ag is the shift of the g- 
factor which, for cuprates, is of order 0.1. Here, however, Vi mp and v so will be treated as 
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independent variables. Note that, since the momenta k and k' are defined in the x-y plane, 
the spin-momentum dependence of the spin-orbit interaction simplifies to: 



[k x k'] • cr al3 = [k x k'] ■ z a 



z 

a(3- 



(4) 



Since the spins have been quantized along the x-axis, the spin-orbit coupling leads to scat- 
tering events always accompanied by spin-flip transitions. 

The following analysis is simplified by introducing the usual four-components field oper- 
ators [0,0: 



c-ki 



i k 



(5) 



From Eqs.(|T]||) it is possible to evaluate the equation of motion of the field operator \l/k in 
the imaginary time r: 

C^k 



dr 



e(k)p 3 ^ k - A(k)p 2 r 2 ^ k + i> 3 r 3 tf k 



i(k-k')-Ri 



k',i 



v im pp3 + i-j^ L [kxk']-ZT 1 



(6) 



where the products piTj are 4x4 matrices acting on the field operators (||). They are 
constructed by treating the Pauli matrices Tj as elements of the Pauli matrices pi as shown 
in the example below: 



/ -tT 2 ^ 



P2T2 



V 



(7) 



ir 2 

Equation (^|) permits to evaluate the equation of motion of the generalized Green's 
function defined as: 



G(k,k';r) = -(T T ^ k (r)¥ kl (0)), 



(8) 



where T T is the r-order operator. It is straightforward to obtain from Eqs. (||[5p the equation 
satisfied by the generalized Green's function in the Matsubara frequencies uj n = (2n + 1)ttT: 



G(k, k'; iu n ) = 5 kik /G (k; iu n ) + G (k; iu n ) 

EE^ (k - k " )R ^(k,k")G(k",k';^„), 



x 

k" i 



(9) 



where 



v(k, k") = v mpP3 + i|[kx k"] • z n, (10) 

and 

G (k; za> n ) = [za;„ - e(k)p 3 - A(k)p 2 r 2 - Ip^-iY 1 , (11) 

is the Green's function in the absence of impurities. 

The average over all the impurity configurations of Eq. (|5|) leads to the averaged Green's 



function G which satisfies the following Dyson equation [10 



G 1 (k,iu n ) = G 1 (k,iuj n ) -E(k,iu n ), (12) 

where S is the electron self-energy resulting from the average procedure. In this paper, v imp 
and v so are assumed to be sufficiently weak to justify a self-consistent Born approximation 
for the self-energy S. Because of the momentum dependence of the spin-orbit interaction, 
the Feynmann diagrams describing the Born approximation do not involve impurity-spin- 
orbit mixed terms and S is given by the diagrams showed in Fig. |l|. Therefore, by using 
Eq . (|l0l) , the self-consistent Born approximation for E reads: 

S(k; iu n ) = n^ V(k, k')G(k', iu n )V(k', k) 

k' 

= n iVi m pJ2p3 G ( k '^ iuJ n)p3 + rii-^Yj\ k X k!\ 2 T X G(k! ] iuJ^Tx, (13) 
k' K F k' 

where is the concentration of impurities. 

Equations ( |12|) and (|H|) must be solved self-consistently and the solution can be written 
in terms of the following renormalized Green's function |7| JT0|JTT|] : 

G _1 (k; iu n ) = i{u - ilp 3 a 3 ) - p 3 (e - ^Ap 3 cr 3 ) - p 2 cr 2 (A - iQp 3 a 3 ), (14) 



where for brevity the momentum and frequency dependence of the tilded quantities has 



been omitted. The renormalized quantities can be calculated by substituting Eq.([L4]) into 
Eqs. (p^JT3| ). If there is particle-hole symmetry, the quasiparticle dispersion remains unaf- 
fected by the presence of impurities, i.e., e = e(k) and A = 0. For the other quantities it is 
useful to introduce the variables u± and A± defined by: 

Cj ± =Cj±U; A± = A±ifi. (15) 
In this way the self-consistent equations become: 



^ ^ |k x k| 2 i^ 
£ V e(k') 2 + A 2 , 

A ± = A(k) + n, Lp E e(kl)a+ % + ^ + 4 E e ^ a X + %% - (17) 

The summations over momenta are transformed in integrations over energy according to the 
usual procedure: 

J (2tt) 2 Jo 2tt J ' K ' 

where N is the electronic density of states per spin state at the Fermi level. Performing the 



integration over the energy e, equations (jig) and ([T7|) reduce to: 

1 /• g> ± 1 rd^ ff(0,0')^ T 
a;+ = c<; n ± 2/ H / ~ — H / = — -, 19 

2tJ 2vr [^ + A±(0') 2 ] 1/2 2tt p| + A T (0') 2 ] 1 / 2 

A ± (0) = A(0) + -I / ^ + JL / ^WlML, ( 20 ) 

MW lW 2ri 2tt p 2 +A ± (0') 2 ] 1 / 2 r so i 2vr pl + A^^O 2 ] 1 ^' 1 ; 

where r -1 and (r^)" 1 are the scattering rates for the non-magnetic and spin-orbit impurities, 

respectively. They are given by: 



1 2 1 

- = 2im i v imp No , - 
r r. 



2imiv1 mp N Q , — = imiV 2 so N . (21) 



In Eqs.( |T9pO| ), the function S((j),<fi') stems from the angular dependence of the spin-orbit 
factor |k x k'| 2 by defining and 0' as the polar angles of the versors k and k', respectively. 
In explicit form, the function S((j), <$') is given by: 



S(<f), 0') = cos(0) 2 sin(0') 2 + sin(0) 2 cos(0') 2 - \ sin(20) sin(20')- (22) 

The presence of such angular function leads to important differences between non-magnetic 
and spin-orbit impurity effects also for zero magnetic field. In fact, non-magnetic impurities 
does not renormalize the gap function when this has ci-wave symmetry || whereas the spin- 
orbit interaction, by means of the angular function S(<f>, <p'), provides a finite renormalization. 
This can be readily seen by realizing that if A(0) is of the form given by Eq.(^), then a 
consistent solution of Eq.(|20|) is provided by setting A±(0) = A±cos(20), where A± is the 
solution of: 

~ . 1 f d(b , ,. 9 AzpCos(2</>) , , 

A± = A + — / — sin (6 f— ±- -, 23 

t so J 2tt W) p| + A 2 f cos(20) 2 ] 1 /2' 1 ; 

and, in the same way, Eq.(|i~9|) becomes: 

1 f d(f) G)± 1 f d(f) . , x2 CJzr 



1 [d<p uj± 1 fd(p . , r 

UJ+ = (Jj„ ± 11 H / = ; — / — sin 

2rJ 2it\lu1 + A! cosf2^ 2 lV2 T sn J 2n KV) 



27r[a;| + A 2 b cos(2(/.) 2 ] 1 /2 r so J 2tt pi + A 2 F cos(20) 2 ] 1 /2 

(24) 



In obtaining Eqs.(E3l24), it has been used the identity: 



/ ^ cos(0) 2 /[cos(20)] = sin(0) 2 /[- cos(20)], (25) 

J ATI J Z7T 

where /[cos(20)] is a general function of cos(20). 

As expected, the scalar impurity scattering contribution has disappeared from the gap 
renormalization (^3|). On the contrary, the spin-orbit interaction modifies the gap function 
because of the presence of the angular function fl22"|) . Moreover, equations ( ^3]) and ([24]) are 
renormalized in a different way by v so so that, even at zero magnetic field, the spin-orbit 
interaction contributes to the thermodynamic and spectral properties of ci-wave supercon- 
ductors. In fact, all the measurable quantities can be expressed in terms of u± = u±/A± 



i~T|| which from Eqs.(^,[24]) satisfies the following equation: 



_ u n ±il 1 fd4 u± 1 f d4_ . 2 u T - u± cos(20) 

U± " A + 2Ar J 2tt [cos(20) 2 + u 2 ± ]W + Ar so J 2n Sm{<P) [cos(20) 2 + u%} l / r 1 ' 



The above equation should be compared with the corresponding expression for the two- 
dimensional isotropic s-wave case which, by setting A(0) = A in Eqs. ( |19| , [20| ) , is found to be 

100 

u n ±il 1 u T - u± 

M± = -E- + 2A^ |TTs|F- (27) 

where the contribution of the impurity scattering has vanished because of Anderson's theo- 



rem. When H = 0, equation (|27|) reduces to w+ = U- — u; ra /A and does not depend on the 
spin-orbit scattering rate while Eq.(^(J) still depends on r and r so . In fact, in a c?-wave super- 
conductor, both the non-magnetic impurity and the spin-orbit scatterings are pair breaking 
and they tend to suppress superconductivity ||. When H ^ 0, it is therefore expected that 
the Zeeman response of a <i-wave superconductor differs qualitatively from that of a s-wave 
condensate. 

III. PHASE DIAGRAM 

Equation (p6|) permits to obtain all the informations needed to calculate the phase dia- 
gram of a dirt (i-wave superconductor in a Zeeman magnetic field. Let us start by considering 
the self-consistent equation for the order parameter A: 



A-^rEE^m^^^^g/gR. ^^ , (28) 



cos(20) 2 



2tt [cos(20) 2 + u\]W 



where Vq is the pairing interaction and A = VoiVo. The summation over the frequencies is 
implicitly assumed to be restricted by a cut-off energy. However, both the cut-off frequency 
and the pairing interaction can be absorbed in the definition of the critical temperature T c o 
for a pure superconductor (r _1 = 0, r" 1 = 0) without external magnetic field. In this way 
the gap equation can be rewritten as: 

In the hypothesis that the transition to the normal state is of the second order (see below), 
the critical temperature T c is obtained from Eq . (|28D by setting A — > and it is given by: 
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4r so fe/ V 2 27rT c/ \ 4r ^/ V 2 2tcT c, 

(30) 



where a = (r 1 + r so 1 )/47rT c and 6 = [l/(4r so ) 2 — J 2 ] 1 / 2 and ■?/> is the di-gamma function. 
When / = [i B H = equation ( j30|) reduces to: 



©=*G)-* {I + 4^> + j2^)' (31) 
which coincides with the result obtained in Ref . || in the weak scattering limit . Equation 



(pl|) shows that, even at zero magnetic field, the spin-orbit scattering contributes together 
with the non-magnetic impurity scattering to the suppression of T c . 

For large enough values of the external magnetic field, the transition to the normal state 
becomes of first order Q. This situation is studied by evaluating the difference of the free 
energy between the superconducting and the normal states AF = F s — F n . If, by rising the 
temperature and/or the magnetic field, AF changes sign while A remains finite, then the 
system undergoes a first order phase transition to the normal state with critical field B c and 



T c determined by AF = 00. By following Ref. ||11|| , AF is obtained by: 

AF= / dV A 2 , (32) 
Jo 

and by using Eqs. (|26"| , [2"8|) one readily finds: 

AF = -iV.A2.rg/gRe + - 2u + - ^^^p } ■ (33) 

The numerical solution of equations ( |30| ) and (^) are shown in Fig. |^ for the pure limit and, 
for comparison, the ci-wave solution is plotted together with the s-wave one. In the phase 
diagram, the solid and dashed lines are solutions of equations fl3Cf ) and (0), respectively. For 
both (i-wave and s-wave, the transition to the normal state is of second order for T/T c q > 0.56 
|1],|2||J while for lower temperatures the transition, marked by dashed lines, becomes of first 
order. For T/T c0 < 0.56 the solid lines represent the supercooling field |T]-§[. As already 
said in the introduction, at zero temperature the first order transition to the normal state is 
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obtained by critical fields [1bH c /A = l/\/2 [[| for s-wave and /i£if c /A — 0.56 for d-wave 
In this paper, the Fulde-Ferrel-Larkin-Ovchinnikov state |nj which appears at low 
temperatures has not been considered since disorder tends to restore the zero momentum 



pairing |H| . For the pure (i-wave case, the reader can find the phase diagram including the 
non-zero momentum pairing state in Ref. [Q. 

Although for the pure limit the phase diagrams of the Zeeman response of s-wave and 
d-w&ve superconductors are qualitatively similar, they drastically differ when the coupling 
to the non-magnetic and spin-orbit impurity scatterings is switched on. In Figs. |] and [|, the 
phase diagrams for s-wave and (i-wave superconductors are plotted for finite values of r _1 and 
r s o . In both figures, the impurity scattering parameter b n = l/(2A r) is set equal to 0.1, 
while the spin-orbit scattering parameter b so = 1/ (2Aor so ) assumes four different values: 
b so = 0, 0.06, 0.12, 0.16. In the s-wave case, Fig. 3, the phase diagram is insensitive to 
b n 7^ 0, while finite values of b so enhance the critical field for all temperatures. Moreover, the 
temperature interval of first-order phase transition (dashed lines) decreases as b so increases 
and for b so > 2.32 the transition becomes continuous for all temperatures ||. This behavior is 
due to the spin-mixing effect of the spin-orbit interaction which lowers the Zeeman response 
and consequently the depairing effect of the magnetic field. On the other hand, in the (i-wave 
case shown in Fig. 4, the spin-orbit scattering is pair breaking and for b so > the critical 
field is lowered. This situation can be understood by realizing that finite values of b so lead to 
a weakening of the superconducting state || so that, with respect to the b so = case, lower 
values of H are needed to suppress completely superconductivity. Another striking feature 
is that, due to the nodes of the (i-wave order parameter, the b so dependence of the phase 
diagram is much stronger than for the s-wave case. In fact, already for b so = 0.16 there is 
not signature for first-order transition whereas for an s-wave superconductor the first-order 
transition disappears only for b so > 2.36, i. e., a difference of one order of magnitude. 

It is important to stress that the remarkable difference between the (i-wave and s-wave 
phase diagrams has been obtained in the Born limit of non-magnetic and spin-orbit impurity 
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scatterings. However it is well known that in high-T c superconductors the effect of disor- 
der is best described by the strongly resonant limit of the impurity potential so that the 
Born approximation may result inadequate. In practice, one should formulate the Zeeman 
response by employing the t-matrix approximation for both the non-magnetic and the spin- 
orbit potentials. Such a calculation has already been reported in Ref. || for zero external 
magnetic field. The generalization for H ^ is currently under investigation and the results 
reported here may provide a useful tool for testing the more general t-matrix solution. 

IV. DENSITY OF STATES 

Having described the effect of b n and b so on the phase diagram, it is interesting to 
investigate also how the spectral properties are modified. To this end, equation ( p6|) must 
be analytically continued to the real axis by setting iu± —>■ u± and iu n — > uj. In this way, 
the quasiparticle density of states (DOS) per spin direction in units of the normal state DOS 
Nq can be calculated by the following expression: 

N±(uj) r d(j) u± 
P±M = AT = sgn w / — Re-^ ,^ 211/2 - 34 

In Fig. |] it is reported the quasiparticle DOS for / = psH = 0.15A , b n = 0.1 and 
different values of the spin-orbit parameter b so . For clarity, the curves with b so ^ have 
been vertically shifted by 0.7, 2 x 0.7 and 3 x 0.7 with respect to those with b so = 0. For 
b so = 0, the two DOS per spin state, p + (dashed lines) and p_ (solid lines), show a clear 
Zeeman splitting and for uj = the total DOS p = p + + p_ is different from zero as expected 
for a c?-wave superconductor. For b so > the total DOS at u = is enhanced at the expenses 
of the coherence peaks which show a decrease of spectral weight. Moreover, at u ~ PbH, 
p_ develops a structure (marked by the arrows) which becomes a peak at b so = 0.16. Such 
a structure is even more visible in Fig. |^ where p± is plotted for b n = 0.1, b so = 0.06 and 
different values of the external magnetic field. The origin of this peak can be understood 
by the following reasoning. At the Fermi wave-vector, and for a pure superconductor, the 
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quasiparticle energies for spin up and down are E±((f>) = A| cos(20)| ± \ibH and therefore, 
depending on the values of H and A, two quasiparticles with different spin orientation 
and angles can have equal energies. For example, for <px — and <p2 = 7f/4, the two 
energies and E + {4>2) are equal to uj = fisH if A = 2/isiJ. Since the spin-orbit 

potential connects quasiparticle states with different spin orientation but equal energies, the 
two states and E + (<p2) are coupled by the spin-orbit interaction and a enhanced 

signal should be expected at u ~ fisH. Note in fact that in Fig. |] the low energy peak is 
more pronounced for \ibH = 0.35A where, since A ~ 0.68A , the condition A = 2/isif is 
nearly fulfilled. 



V. CONCLUSIONS 

In conclusion, it has been shown within the Born approximation that the presence of 
impurity and spin-orbit scattering centres strongly affects the Zeeman response of a d- 
wave superconductor. Both the phase diagram and the quasiparticle density of states show 
features qualitatively different from those of a s-wave superconductor. In fact, by increasing 
the value of the spin-orbit scattering parameter b so = l/(2A r so ) the critical field H c is 
strongly lowered whereas in a s-wave superconductor H c increases. Moreover, the influence of 
b so on the superconducting state is much stronger for the ci-wave symmetry. Concerning the 
spectral properties, the Zeeman splitted density of states of a <i-wave superconductor shows 
interesting features which are missing in a s-wave superconductor. In fact, for sufficiently 
large values of b so and/or H a resonant peak develops at energies close to fisH. The origin 
of this feature is given by the anisotropy of the order parameter and the spin-flip transitions 
due to the spin-orbit scattering. An important open question concerns with the possibility 
of going beyond the Born approximation and employing a t-matrix approach for the Zeeman 
response in the presence of impurity and spin-orbit scattering centres in order to test the 
solidity of the results here presented. 
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FIG. 1. Feynmann diagrams for the self-energy in the self-consistent Born approximation. The 
impurity and spin-orbit interactions are represented by dashed and dot-dashed lines, respectively. 
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FIG. 2. Phase diagram for pure s-wave and ci-wave superconductors in the presence of a Zeeman 
magnetic field. Ao and T c q are the order parameter and the critical temperature without the 
external magnetic field, respectively. For T/T c q > 0.56 the solid lines are the second-order phase 
boundary between the normal (above the solid lines) and the superconducting (below the solid 
lines) states. For T/T c q < 0.56 both the s-wave and the d-wave states show a first-order transition 
to the normal state marked by the dashed lines. In this region, the solid lines represent the 
supercooling fields. 
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FIG. 3. Phase diagram for a s-wave superconductor with impurity and spin-orbit scattering 
centres. The critical field is unaffected by the impurity potential while it increases by increasing 
the spin-orbit scattering parameter b so = 1/(2Aot so ) where Ao is the order parameter in the pure 
limit without magnetic field. The solid and dashed lines have the same meaning as in Fig.|2|. 
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FIG. 4. Phase diagram for a d-w&ve superconductor with impurity and spin-orbit scattering 
centres. The impurity scattering parameter is b n = 1/(2Aot) =0.1 where Ao is the order parameter 
in the pure limit without magnetic field. The solid and dashed lines have the same meaning as in 
Fig|2|. Note that, contrary to the case shown in Fig.||, already for b so = 0.16 the transition to the 
normal state is of second-order for all the temperature range. 
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FIG. 5. Zeeman splitted quasiparticle density of states p+{ui) (dashed lines) and P-{uj) (solid 
lines) for a ci-wave superconductor with b n = 0.1, ^bH/Aq = 0.15 and different values of the 
spin-orbit scattering parameter b so . The curves for different values of b so are vertically shifted by 
multiples of 0.7. Note the structure (marked by the arrows) at uj ~ pbH which develops as b so 
increases 
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FIG. 6. Zeeman splitted quasiparticle density of states p+{u) (dashed lines) and P-(ui) (solid 
lines) for a <i-wave superconductor with b n = 0.1, b so = 0.06 and different values of the external 
magnetic field. The curves for different values of H are vertically shifted by multiples of 0.7. The 
arrows indicate the resonant structure which develops a peak for psH > 0.3 (see text). 
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